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O-i Abstract 

This paper reports the evaluation of the tenth-order QED contribution to the lepton g — 2 from 

, the gauge-invariant set of 2072 Feynman diagrams, called Set IV, which are obtained by inserting a 

>' 

\Qi second-order lepton vacuum-polarization loop into 518 eighth-order vertex diagrams of four-photon 

<n : 

exchange type. The numerical evaluation is carried out by the adaptive-iterative Monte-Carlo 

(N ; 

. integration routine VEGAS using the FORTRAN codes written by the automatic code-generating 

t-h ; 

t-h I algorithm gencodeA/. Some of the numerical results are confirmed by comparison with the values 

of corresponding integrals that have been obtained previously by a different method. The result 

•rH . 

^ ■ for the mass-independent contribution of the Set IV to the electron g — 2 is —7.7296 (48)(a/7r) 5 . 

b : 

There is also a small mass-dependent contribution to the electron g — 2 due to the muon loop: 
-0.01136 (7)(a/vr) 5 . The contribution of the tau-lepton loop is -0.0000937 (104) (a/vr) 5 . The 
sum of all these contributions to the electron g— 2 is —7.7407 (49)(a/-7r) 5 . The same set of diagrams 
enables us to evaluate the contributions to the muon g — 2 from the electron loop, muon loop, and 
tau-lepton loop. They add up to —46.95 (17)(a/7r) 5 . 

PACS numbers: 
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I. INTRODUCTION 



The anomalous magnetic moment g — 2 of the electron has played the central role in 
testing the validity of quantum electrodynamics (QED) as well as the standard model. On 
the experimental side, the latest measurement of a e = (g — 2)/2 by the Harvard group has 
reached the precision of 0.24 x 10~ 9 jl, 2|: 



a e (HV08) = 1 159 652 180.73 (0.28) x 10~ 12 [0.24ppb] . 



The theoretical prediction thus far consists o 
5j, direct evaluation of hadronic corrections 



QED corrections of up to the eighth order 



13 



(1) 

EL 



and electroweak corrections scaled 
15|. To compare the theory with the 



down from their contributions to the muon g — 2 
measurement (JT]), we also need the value of the fine structure constant a determined by 
a method independent of g — 2 . The best value of such an a available at present is one 
obtained from the measurement of h/m-Rb, the ratio of the Planck constant and the mass of 
Rb atom, combined with the very precisely known Rydberg constant and mRb/m e : [16| 



a -1 (RblO) = 137.035 999 037 (91) [0.66ppb]. 



With this a the theoretical prediction of a e becomes 



a e (theory) = 1 159 652 181.13 (0.11)(0.37)(0.02)(0.77) x 10 



(2) 



-12 



(3) 



where the first, second, third, and fourth uncertainties come from the calculated eighth-order 
QED ten* fl, a crude tenth-order estimate Q , the hadronie and e.ectroweah contributions, 
and the fine structure constant (J2J), respectively. The theory ([3]) is in good agreement with 
the experiment ([T|): 



a e (HV08) - a e (theory) = -0.40 (0.88) x 10 



-12 



(4) 



proving that QED (standard model) is in good shape even at this very high precision. 

An alternative and more sensitive test of QED is to calculate a from the experiment and 
theory of g— 2 , both of which have very high precision, and compare it with a _1 (Rbl0). 
The experiment and theory of the electron g — 2 leads 



a -1 (a e 08) = 137.035 999 085 (12)(37)(2)(33) [0.37ppb], 



(5) 
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where the first, second, third, and fourth uncertainties come from the eighth-order QED 
term, the tenth-order estimate, the hadronic and electroweak contributions, and the mea- 
surement of a e (HV08), respectively. 

Although the uncertainty of a _1 (a e 08) in (jSJ) is almost a factor 2 smaller than that of 
a _1 (RblO), it is not a firm factor since it depends on the estimate of the tenth-order term, 
which is only a crude guess [17|. For a more stringent test of QED, it is obviously necessary 
to evaluate the actual value of the tenth-order term. To meet this challenge we launched 



several years ago a systematic program to evaluate the complete tenth-order term [18H20]. 

The tenth-order QED contribution to the anomalous magnetic moment of an electron 
can be written as 



1 7T 



Af 0) + A^°\m e /m^) + A ( 2°\m e /m T ) + A^°\m e /m^m e /m T ) , (6) 



where the electron-muon mass ratio m e /m M is 4.836 331 66 (12) x 10 -3 and the electron-tau 
mass ratio m e /m T is 2.875 64 (47) x 10 [121]. In the rest of this article the factor (a/ V) 5 
will be suppressed for simplicity. 

The contribution to the mass-independent term can be classified into six gauge- 

invariant sets, further divided into 32 gauge- invariant subsets depending on the nature of 
closed lepton loop subdiagrams. Thus far, the numerical results of 29 gauge-invariant sub- 



sets, which consist of 3856 vertex diagrams, have been published [3|, l2l|-|27l]. Five of these 



28. 



i\) gaug 
QBE 



29). They are in good agreement with our 



29 subsets were also known analytically 
calculations. 

In this paper we report the result of evaluation of A^ from the set, called Set IV, 
which consists of 2072 Feynman diagrams. Sec. [Til outlines our formulation of Feynman- 
parametric integrals of Set IV. Sec. II III presents the residual renormalization formula, which 
summarizes the result of derivation described in detail in Appendix [A] Numerical results 
for several cases of mass dependence are described in Sees. HVllVl and[VTl Sec. IVHI discusses 
the results obtained in this paper. 

II. CONSTRUCTION OF FEYNMAN-PARAMETRIC INTEGRALS 

All 2072 diagrams of Set IV can be derived from the 518 eighth-order diagrams of four- 



photon-exchange type 



30] , called Group V, by inserting a second-order vacuum-polarization 
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FIG. 1: The eighth-order Group V diagrams. The solid line represents the electron propagating 
in a weak constant magnetic field. 

loop in the photon lines of Group V diagrams in all possible ways. In practice, we have 
therefore to deal with only 518 diagrams. This can be reduced further to the 47 self-energy- 
like diagrams of Fig. [T] as follows. 

Let A u be the sum of 7 vertex diagrams that are obtained from any self-energy-like 
diagram of Fig. [I] by inserting a magnetic vertex 7 1 ' in all possible ways. The set 
of these vertex diagrams, taking account of doubling due to time-reversal, represents the 
original 518 vertex diagrams. The next step is to rewrite this A u as 



A"(p, q) ^ -q" 



dA^(p,q) 
dq„ 



g=0 



0£(p) 



(7) 



for small q, with the help of the Ward-Takahashi(WT) identity, where p — q/2 and p + q/2 
are the 4-momenta of incoming and outgoing lepton lines and (p — q/2) 2 = (p + q/2) 2 = m 2 . 
The g — 2 term is projected out from the right-hand side of Eq. ([7]). 

The properties of the Feynman-parametric integrals corresponding to the diagrams of 
Fig. [T] have been studied and described in detail in Ref. |5j. Each diagram Q of Fig. [T] is 
represented by a momentum integral using the Feynman-Dyson rule. Introducing Feynman 
parameters z%, Z2, ■ ■ ■ , for the electron propagators and z a , z&, z c , Zd for the photon prop- 
agators, we carry out the momentum integration analytically using a home-made program 



written in FORM 3l|, which gives an integral of the form 



M c 



(dz) 



1 / Eq + Co Ei + G\ 



U 2 V 3 



U 3 V 2 



+ ■ 



+ 



N + Z Ni + Zi _ 



U 2 V 



U 3 V 3 



(8) 



where E n , C n , N n and Z n are functions of Feynman parameters. The subscript n of E n , etc., 
indicates that it is the n contraction terms of diagonalized loop momenta and proportional 
to the product of n factors of B^s. The "symbolic" building blocks Ai, B^, C^, for i,j = 
1,2, ... ,7 are also functions of Feynman parameters. U is the Jacobian of transformation 
from the momentum space variables to Feynman parameters. V is obtained by combining 
denominators of all propagators into one with the help of Feynman parameters. It has a 
form common to all diagrams of Fig. [TJ 

7 d 

v = 5>(i-^K + ]T^, ( 9 ) 

i=l k=a 

where and are the rest masses of electron i and photon k, respectively. Ai is the scalar 
current defined by 



A; 



and 



(dz)g = Y\ dziS(l - ^ Zj). 



(10) 



(11) 



ieg 



See, for example, Ref. |32( for definitions of B^ and Cy. The form of Ai as a function of 
Feynman parameters depends on the structure of individual diagram. However, as is shown 
in Eq. (Q, the expression of V in terms of Ai is identical for all diagrams of Fig.[H Individual 
diagram of Fig. [1] will be denoted as Mg and their assembly will be collectively denoted as 
M 8 . 

We have developed two independent sets of numerical programs of Mg based on the WT- 
summed amplitudes. The first formulation was developed in 1970's and given in Ref. [33 ]. 



The second formulation used the automation code GENCODE./V [19j, [20j. The unrenormalized 
amplitudes and the UV-subtraction terms are the same, but the IR-subtractions are slightly 
different in two formulations. The detail of UV- and IR-subtraction terms in the second 
formulation is briefly described in Sec. Ill IL After taking account of the difference in two 
formulations, the equivalence of two formulations is established 20] . Once we have the cor- 



rect programs of the eighth-order Group V diagrams, the insertion of a vacuum-polarization 



M 47 M 4ZP2 

FIG. 2: The eighth-order diagram M47 of Group V and the tenth-order diagram M^j t P2 of Set 
IV. The diagram M47 P2 represents the sum of diagrams obtained by inserting a second-order 
vacuum-polarization loop into each of four photon lines of the eighth-order diagram M47. 

loop is an easy task to carry out. Fig. [2] shows a typical self-energy- like diagram of the 
tenth-order Set IV. 

As is well-known, the insertion of a vacuum-polarization loop in an internal photon line 
can be expressed as a superposition of massive vector particle propagators. In other words 
all we have to do is to replace the mass square A 2 of one of the photons in Eq. ([9]) by pit): 

\ 2 ^p(t) = ^ (12) 

where m vp is the mass of the fermion forming the vacuum-polarization loop, to multiply the 
resulting eighth-order integral with the spectral function 

"sW = H< 2 ) ' (13) 

and to integrate over the interval < t < 1. 



This is easy to implement in the second formulation |19l . l20j ] since the function V is 
unambiguously identifiable. Unfortunately, in the first formulation [33], it is difficult to 
implement this procedure for some diagrams because the " denominator function V" was 
used to replace parts of numerators in order to reduce the size of integrands and accelerate 
the computing speed. For this reason, it is difficult to apply Eqs. (fT2"j) and (fl3"|) to these 
integrals. Thus, direct comparison of two methods is feasible only for those of Set IV 
diagrams in which the function V can be clearly distinguished from other terms of the 
numerator. We therefore report here only the results of the second formulation. Since 

n 

the equivalence of two methods has been well established [5|], this does not diminish the 
reliability of our numerical results. 
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III. RESIDUAL RENORMALIZATION 



In our approach based on numerical integration the integrals of individual diagrams must 
be made convergent before they are integrated numerically. This is achieved in the following 
manner. 

Suppose the integral Mg has a UV divergence arising from a subdiagram S. Then we 
construct another integral K$Mg by applying a i^-operation, which identifies and extracts 
the UV divergent part of Mg by a simple power counting rule. This integral has the following 
properties: 

• It has the same domain of integration and the same UV divergence as Mg. Thus it 
subtracts the UV divergence of the latter point-by-point in the domain of integration. 

• If S is a vertex diagram, i^-operation Kg on Mg factorizes exactly into the product of 
lower-order quantities as 

K S Mg = L^Mg /S . (14) 

If S is a self-energy diagram, if-operation K5 on Mg turns exactly into the sum of two 
terms of the form 

K s M g = SmfMo/sv) + BfM g/s( , y (15) 

Here L^ v , B^ y , and drn^ are UV-divergent parts of the renormalization constants 
Ls, Bs, and 8ms- Mg/s is the magnetic moment corresponding to the diagram Q/S 



32] for further 



obtained by shrinking the subdiagram S of Q to a point. See Ref. 
details. 

An IR divergence of Mg arises from a subdiagram T that is the reduced diagram T = Q/S 
of a self-energy subdiagram S of the diagram Q. In this case we run into two kinds of IR 
divergence. One arises when a self-energy subdiagram S behaves as a self-mass term. The 
standard mass-renormalization on Q subtracts 5msMg/ s ^*^ from Mg while Kg operation 
of Eq. (fT5"|) subtracts 5m^ Y Mg/su*)- Thus, after subtraction by the K$ operation, we are 
left with (5ms — 5m s rv )Mg/s(i*), which has a linear IR divergence because of divergent 
Mg/su*), except when 8ms = <5m^ v . The easiest way to deal with this problem is to subtract 
8ms entirely instead of only 5m^ v . We call this .^-subtraction, which is incorporated in 

GENCODEiV. 
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The other IR divergence occurs when a self-energy-like subdiagram S behaves as a mag- 
netic moment amplitude. The remaining diagram T can be mimicked by a vertex diagram by 
shrinking the subdiagram S to a point. This divergence is only logarithmic and the subtrac- 
tion term can be constructed by applying the /-subtraction I7- on the UV-finite amplitude 



Mg, which is shown to factorize as 



20] 



l T M g = L*M 5 , (16) 

where is the part of the vertex renormalization constant Lj- that remains after all UV- 
divergent pieces are subtracted out. 

These operations, carried out for all divergent subdiagrams of the unrenormalized integral 
Mg, create a UV-finite and IR-finite integral AMg. For a full account of these operations 



see Refs. 



19 



20]. 



Since this scheme is different from the standard on-the-mass-shell renormalization, it is 
necessary to make an adjustment, called residual renormalization, which accounts for the 
difference of the standard renormalization and the UV-divergent (and IR-divergent) parts 
generated by if-operation (and i/iZ-subtractions). 

The residual renormalization terms of individual diagrams must then be summed up 
over all diagrams involved. As the order of perturbation increases the total number of 
terms contributing to the residual renormalization increases rapidly so that it will become 
harder and harder to manage. Fortunately, the sum of all residual terms can be expressed 
concisely in terms of magnetic moments and finite parts of renormalization constants of 
lower orders, whose structure is closely related to that of the standard on-the-mass-shell 
renormalization. This observation enables us to obtain the sum of residual renormalization 
terms of all integrals starting from the expression of the standard renormalization. This 
approach is described in detail in Appendix |A] for the eighth-order g — 2 after simpler cases 
of fourth- and sixth-orders are described for illustration of our method. 

Since diagrams of Set IV are obtained from the magnetic moment contribution M 8 of 
518 eighth-order vertices of four- photon-exchange type by inserting a second-order vacuum- 
polarization subdiagram in all possible ways, the residual renormalization term of the Set IV 
is readily derived from that of the residual renormalization term of Mg. Namely, insertion 
of a closed loop of the lepton l 2 in the internal photon lines of Group V diagrams of lepton 
l\ given in Eq. (IA35j) in all possible ways leads to the renormalized contribution of Set IV 
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to the lepton g — 2 of the form: 

A (10) [Set IV (Wa) ] = AM^pz 

- hAM^ALB 2 

- 5AM 6 ALB^ h pf 

+ AAfiJ^ ) (-3ALB 4 + 9(ALB 2 ) 2 ) 

+ AM 4 (-3ALB^ + 18ALB 2 ALB^) 

+ Ml P h 2 \-ALB 6 + 6ALB4ALB2 - 5(ALB 2 ) 3 ) 

+ M 2 (-ALB^ + 6ALB^p2 ALB 2 

+ 6ALB A ALB%$ - lh{ALB 2 f ALB { £$) 
+ M 2 ( ^ ) A5m 4 (4AL 2 * + AB 2 *) 
+ M 2 A8m1^p2 (4AL 2 * + AB 2 ») 

+ M 2 A5m 4 (4AL^|, ) 2 + AB$ h p \), (17) 

where superscripts such as and (/ 2 / 2 ) are omitted for terms which are independent of 
rest mass. See Appendix |A] for the explanation of notations. 

AM^p l 2 is the sum of 74 WT-summed integrals enhanced by the insertion of vacuum- 
polarization- loop. Each of these 74 integrals is finite by our construction. Individual terms of 
residual renormalization are also UV- and IR-finite by construction. Eq. fflTl) thus maintains 
that /l( 10 )[Set IV^ 1 ' 2 - 1 ], which represents the quantity renormalized in the standard manner, 
can be expressed as the sum of completely finite quantities, each of which can thus be 
integrated by numerical means. 

We should like to emphasize that Eq. ffTTI) is analytically exact and involves no approxi- 
mation as far as the subtraction term factorizes exactly as in Eqs. f|T4l) . f|T5|) . and ffT6l) . 



IV. NUMERICAL EVALUATION OF a[ 10) [Set IV] 

AM a p 2 , which is made UV-finite by ff-operation and IR-finite by //^-subtractions, is 



integrated numerically by the adaptive Monte-Carlo integration routine VEGAS [34|. The 
result for (7i/ 2 ) = (ee) are listed in Tables HI and HT1 Auxiliary quantities needed for carrying 
out the residual renormalization are listed in Table [TTT1 Notations are those of Eq. ( fl7l) . 
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TABLE I: Contributions of diagrams Mqi^,. .., -^24^2 of Set IV to a e for (/1/2) = (ee). The 
multiplicity rip is the number of vertex diagrams represented by the integral and is incorporated in 
the numerical value. First 50 iterations are carried out using 1 x 10 8 sampling points per iteration. 
The integrations are continued with 1 x 10 9 sampling points per iteration and iterated as given in 
the second number of the fifth column. The integrals Mi2,p2, -^i6,P2> and Mx^,P2 are evaluated 
with quadruple precision. All other integrals are evaluated with double precision. 



Integral 


Up 


Value (Error) 


Sampling per 


No. of 






including np 


iteration 


iterations 



Mdi,P2 


28 


-0.509 62 


(38) 


1 x 10 8 , 


1 x 


10 9 


50, 


180 


Mq2,P2 


56 


0.060 41 


(98) 


1 x 10 8 , 


1 x 


10 9 


50, 


200 


M)3,P2 


28 


0.829 28 


(60) 


1 x 10 8 , 


1 x 


10 9 


50, 


200 


Mqa,P2 


56 


1.497 37 


(126) 


1 x 10 8 , 


1 x 


10 9 


50, 


360 


Mq5,P2 


56 


0.130 36 


(46) 


1 x 10 8 , 


1 x 


10 9 


50, 


180 


M)6,P2 


56 


-1.084 60 


(94) 


1 x 10 8 , 


1 X 


10 9 


50, 


220 


Mo7,P2 


56 


-1.178 02 


(48) 


1 x 10 8 , 


1 X 


10 9 


50, 


180 


Mos,P2 


56 


-1.415 81 


(121) 


1 x 10 8 , 


1 X 


10 9 


50, 


360 


Mo9,P2 


56 


-0.011 71 


(95) 


1 x 10 8 , 


1 X 


10 9 


50, 


360 


Afio,P2 


56 


-0.816 12 


(107) 


1 x 10 8 , 


1 X 


10 9 


50, 


360 




Zo 


U. 1 05 1 O 


(^o ) 


1 X 1U , 


1 v 

1 X 


in 9 


ou, 


zuu 


M\2,P2 


28 


-1.631 37 


(37) 


1 x 10 8 , 


1 X 


10 9 


50, 


5 


A^13,P2 


56 


-2.353 59 


(45) 


1 x 10 8 , 


1 X 


10 9 


50, 


200 


Ml4 ; P2 


56 


0.685 64 


(83) 


1 x 10 8 , 


1 X 


10 9 


50, 


200 


M 15i p 2 


56 


0.461 55 


(43) 


1 x 10 8 , 


1 X 


10 9 


50, 


200 


A^16,P2 


56 


1.763 95 


(111) 


1 x 10 8 , 


1 X 


10 9 


50, 


025 


Mi7 ; P2 


56 


3.290 90 


(120) 


1 x 10 8 , 


1 X 


10 9 


50, 


340 


A^18,P2 


56 


-0.052 73 


(53) 


1 x 10 8 , 


1 X 


10 9 


50, 


5 


Afig ; P2 


28 


-1.403 73 


(6) 


1 x 10 8 , 


1 X 


10 9 


50, 


180 


A^20,P2 


56 


0.856 32 


(43) 


1 x 10 8 , 


1 X 


10 9 


50, 


180 


A^21,P2 


28 


0.360 89 


(5) 


1 x 10 8 , 


1 X 


10 9 


50, 


180 


A^22,P2 


56 


-0.743 60 


(46) 


1 x 10 8 , 


1 X 


10 9 


50, 


180 


M 2 3,P2 


56 


-1.120 08 


(90) 


1 x 10 8 , 


1 X 


10 9 


50, 


200 


M 2 4,P2 


56 


0.870 63 


(59) 


1 x 10 8 , 


1 X 


10 9 


50, 


200 



Substituting these quantities in Eq. (117]) we obtain 

4 10) [Set IV (ee) ] = -7.729 6 (48) 



(18) 



V. NUMERICAL EVALUATION OF A\ w> (m e /m M ) AND A\ u \m e /m T ) 

Once FORTRAN programs for mass- independent contributions are obtained, it is straight- 
forward to evaluate the contribution of mass- dependent terms such as A^ \m e /m fl ). We 
simply have to choose an appropriate rest mass for the loop fermion l 2 . The result for 
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TABLE II: Contributions of diagrams M25 P2,- ■ ■ , M47 p 2 of Set IV to a e for (ZiZ 2 ) = (ee). The 
multiplicity rep is the number of vertex diagrams represented by the integral and is incorporated 
in the numerical value. All integrals are evaluated with double precision. 



Integral 


np 


Value (Error) 
including np 




Sampling per 
iteration 


No. of 
iterations 


-^25,P2 


28 


-0.696 


52 


(25) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


180 


-^26,P2 


28 


-0.432 


10 


(51) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


200 


M 2 7,P2 


56 


1.120 


35 


(87) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


200 


A^28,P2 


56 


0.783 


12 


(94) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


240 


M29,P2 


28 


1.495 


00 


(92) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


200 


M^P2 


28 


-0.850 


74 


(95) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


200 


A^31,P2 


28 


2.297 


81 


(13) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


180 


A^32,P2 


56 


-2.675 


78 


(35) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


180 


A^33,P2 


28 


-0.960 


21 


(6) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


180 


A^34,P2 


56 


-0.967 


04 


(34) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


180 


-^35^2 


56 


-0.796 


99 


(36) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


180 


A^36,P2 


56 


1.171 


39 


(41) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


180 


M37 i P2 


28 


0.709 


94 


(13) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


180 


A^38,P2 


28 


0.247 


72 


(29) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


200 


Af39,P2 


56 


-0.830 


00 


(30) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


180 


M 40i P2 


56 


-0.499 


07 


(47) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


200 


M4i ; P2 


28 


-1.083 


44 


(71) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


200 


M 42i P2 


28 


0.576 


12 


(76) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


200 


M 43i p 2 


28 


-1.074 


51 


(41) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


200 


M44 ; p 2 


56 


1.919 


47 


(60) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


200 


-^45^2 


28 


0.011 


51 


(37) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


200 


M46,P2 


28 


-0.588 


88 


(73) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


200 


A^47,P2 


28 


-0.102 


58 


(65) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


200 



A 2 \m e /m^) is listed in Tables HVl and IVl From these Tables and the additional data listed 
in Table |Vl] we obtain 

4 10) [Set IV (em) ] = -0.011 36 (7). (19) 

We have also computed the contribution of tau-particle loop A^2°\m e /m T ) } which we give 
without details: 

4 10) [Set IV (et) ] = -0.000 093 7 (104). (20) 

The contribution of the muon loop ffT9l) is about 0.13 % of the electron loop contribution 
( TT8]) . while the contribution of the tau-lepton loop (|20|) is much smaller than the uncertainty 
of ffTHj) and hence completely negligible at present. 
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A T3T T? TTT. 

IAdLHj 111: 

Notations are 


Residual renormalization constants 
tnose 01 H/q. i|i < p. 


needed for the calculation 


ot de petlv v y J. 


AM 6jP2 


1.014 060 (30) 


AM 6 


0.425 820 (14) 


AM 4 , P2 


-0.106 707- •• 


AM 4 


0.030 833 •• • 


M 2 ,p 2 


0.015 687- • • 


M 2 


0.5 


ALB S} p 2 


0.351 54 (93) 


ALB 6 


0.100 86 (77) 


ALB^P 2 


-0.114 228 (17) 


ALB 4 


0.027 930 (27) 




u.d /» i oy [Lo ) 


A(5?7J4 


i ona Q/in ( r >i\ 


ALB 2j p 2 


0.063 399 •• • 


ALB 2 


0.75 


AZy2* ) P 2 


-0.023 531 •• ■ 


AL 2 * 


-0.75 


Ai?2*,P2 


0.047 062 •• • 


AB 2 * 


1.5 



VI. CONTRIBUTION TO THE MUON g-2 

The muon g—2 also receives contributions from the Set IV. The contributions coming from 
the electron loop (l±l 2 ) = (me) are listed in Tables IVTll and rVTTTl Auxiliary quantities needed 
to carry out the residual renormalization are listed in Table I VI I From these quantities we 
obtain 

4 10) [Set IV (me) ] = -38.79 (17). (21) 

We also obtained the contribution of the tau-lepton loop A% (m e /m T ). The result is listed 
in Tables HXl and IXl From these Tables and the additional data listed in Table [VTI we obtain 

4 10) [Set IV (m4) ] = -0.435 7 (25). (22) 

Including the mass-independent contribution (fTBl . the total contribution to the muon g — 2 
amounts to 

A^Pet iy(™e+mm+mt)^ = _ 46 _g 5 (23) 

VII. DISCUSSION 

Since the reliability of the eighth-order term M§ is crucial for the validity of our work on 
the Set IV, let us sketch briefly how we established the validity of M 8 . See Ref. [5j for detailed 
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TABLE IV: Contributions of diagrams Moi P2,- • • , M24P2 of Set IV to a e for (J1I2) = (em). The 
multiplicity rip is the number of vertex diagrams represented by the integral and is incorporated in 
the numerical value. The integral M\2 t p2 is evaluated with quadruple precision. All other integrals 
are evaluated with double precision. 



Integral np 



(em) 



Value (Error) 



Sampling per 



No. of 





including np 




iteration 




iterations 


28 


0.000 


759 


(3) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


56 


0.000 


205 


(8) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


28 


0.001 


757 


(6) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


56 


0.001 


170 


(15) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


56 


0.000 


197 


(7) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


56 


-0.001 


845 


(12) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


56 


-0.001 


866 


(5) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


56 


0.000 


922 


(11) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


56 


-0.001 


077 


(16) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


56 


-0.001 


316 


(9) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


28 


0.000 


310 


(3) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


28 


0.000 


822 


(1) 

\ J 




1 X 


10 7 




20 


56 


-0.001 


434 


(8) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


56 


0.000 


301 


(10) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


56 


0.000 


141 


(5) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


56 


0.000 


920 


(8) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


56 


0.002 


263 


(14) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


56 


0.000 


285 


(4) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


28 


-0.001 


671 


(1) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


56 


0.000 


738 


(12) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


28 


0.000 


210 


(1) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


56 


0.000 


880 


(8) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


56 


0.000 


105 


(27) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 


56 


0.000 


371 


(11) 


1 


x 10 7 , 


1 X 


10 8 


100, 20 



M, 



01,P2 



M, 



(em) 



02,P2 



M, 



(em) 



03,P2 
M (em) 
iw 04,P2 
M (em) 

M (em) 

JW 06,P2 

M (em) 

M (em) 

JW 08,P2 

M (em) 

JW 09,P2 

M (em) 

IV1 W,P2 

M (em) 

M (em) 
1V1 12,P2 
M (em) 
iw 13,P2 
(em) 
1V1 IA,P2 
M (em) 

M (em) 

M (em) 

M (em) 
JW 18,P2 
(em) 
IV1 19,P2 
r (em) 
20,P2 



M 



M (em) 
JW 21,P2 

mH 

JW 22,P2 
M (em) 
iw 23,P2 
r (em) 
24,P2 



M 



accounts. Our approach was to evaluate the diagrams contributing to M s in two independent 
ways. The first method is to apply the scheme formulated more than 30 years ago 33[. The 
revised numerical evaluation by this formulation was reported recently 0, 0]. The second 
approach relies on the FORTRAN codes written by the automatic code-generator gencodeV 



[19. 



201 ] . This method treats the self- mass renormalization terms and IR divergent terms 



differently from the first method so that they can be regarded as practically independent 
of each other. Comparison of the results of these two methods revealed that the first one 
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TABLE V: Contributions of diag: 
multiplicity np is the number of 
in the numerical value. All integr 



rams Mi^pi r . . , M^ } P2 of Set IV to a e for (hfo) = (em). The 
vertex diagrams represented by the integral and is incorporated 
als are evaluated with double precision. 



Integral n p 



(em) 



Value (Error) 
including np 



Sampling per 
iteration 



_ No. of 
iterations 



M, 



25,P2 



M, 



(em) 



26,P2 
M (em) 
1V1 27,P2 
M (em) 
IV1 2S,P2 
M (em) 
iw 29,P2 
M (em) 

M (em) 

M (em) 

1V1 32,P2 

M (em) 

IV1 33,P2 

M (em) 

1V1 3A,P2 

M (em) 

M (em) 

M (em) 

M (em) 

M (em) 

IV1 39,P2 

M (em) 

JW 40,P2 

M (em) 

JW 41,P2 

M (em) 

lvl i2,P2 

M (em) 

iw 43,P2 

M (em) 

M (em) 

JW 45,P2 

M (em) 

iw 46,P2 

M (em) 

iw 47,P2 



98 


— n nm 

u. uu X 


08Q 




i 

X 


v 


10 7 

XU , 


1 

X 


V 


in 8 

xu 


inn 

xuu ■ 


9n 


98 


n nnn 


947 

ill 1 




i 

X 


v 


in 7 

1U . 


1 

X 


v 


10 8 

xu 


1 nn 

xuu ; 


9n 

iU 


56 


n nm 

u.uu X 


937 


f1 51 


i 

X 


ss 


in 7 

xu . 


1 

X 




10 8 

xu 


1 nn 

xuu ; 


9f) 

iU 


56 


n nnn 


nnn 

uuu 




1 

X 


ss 


in 7 

xu . 


1 

X 


v 


10 8 

xu 


1 nn 

xuu ; 




98 


n nnn 

U . UUU 


786 


flM 


1 




10 7 

1 


1 


\<; 


10 s 

X u 


1 nn 


9n 

ZjU 


98 


n nnn 


m a 


(A) 


i 

X 


V 


m 7 

xu . 


1 

X 


S/ 


10 s 

xu 


1 nn 

xuu ; 


91) 


98 


n nm 

u.uu X 


nss 

uoo 




i 

X 


V 


in 7 

xu . 


1 

X 


X 


10 8 

xu 


1 nn 

xuu ; 


9n 

iU 


56 


— n nn2 


282 




I 




10 7 

XU • 


1 




10 8 

xu 


1 nn 

1UU) 


9f) 

iU 


98 


n nnn 

u. uuu 


774 


fl") 


i 

X 




in 7 

xu . 


1 

X 




10 8 

xu 


1 nn 

xuu ; 


9f) 

iU 


UU 


—0 nm 

u.uux 


448 


I'J 


1 

_L 


V 


m 7 

Xu , 


1 

X 


v 


in 8 

xu 


1 nn 

xuu, 






o nnn 

u. uuu 


2Q6 

— • ■ / U 




1 




10 7 

XU • 


1 




in 8 

xu 


mn 

1UU) 




56 


n nm 

u. uu X 


941 


l°J 


i 

X 


V 


in 7 

xu . 


1 

X 




in 8 

xu 


1 nn 

xuu ■ 


91) 


98 


n nnn 

u. uuu 


491 


v°J 


i 

X 


ss 


in 7 

xu . 


1 

X 


N/ 


in 8 

xu 


1 nn 

xuu ; 


9f) 


zo 


n nnn 

u.uuu 


uuO 




1 
1 




m 7 


1 

1 


A 


in 8 

xU 


1 nn 


90 
zu 


56 


0.000 


892 


(5) 


1 


X 


io 7 , 


1 


X 


10 8 


100, 


20 


56 


0.000 


409 


(3) 


1 


X 


10 7 , 


1 


X 


10 8 


100, 


20 


28 


0.000 


775 


(7) 


1 


X 


10 7 , 


1 


X 


10 8 


100, 


20 


28 


0.000 


079 


(4) 


1 


X 


10 7 , 


1 


X 


10 8 


100, 


20 


28 


0.000 


246 


(11) 


1 


X 


io 7 , 


1 


X 


10 8 


100, 


20 


56 


0.000 


546 


(10) 


1 


X 


io 7 , 


1 


X 


10 8 


100, 


20 


28 


0.000 


117 


(2) 


1 


X 


io 7 , 


1 


X 


10 8 


100, 


20 


28 


0.000 


674 


(8) 


1 


X 


10 7 , 


1 


X 


10 8 


100, 


20 


28 


0.000 


203 


(4) 


1 


X 


10 7 , 


1 


X 


10 8 


100, 


20 



had a subtle inconsistency in the handling of some IR subtraction terms. Correcting this 
error we now have two independent evaluations of Mg which agree with each other within 
the precision of numerical integration [5]. 

Although we have not shown the analytic equivalence of the two methods directly, we are 
fully convinced that they are indeed equivalent by proving that they agree to 13 or 14 digits 
(in double precision) at all arbitrarily chosen points in the domain of integration. Only last 
few digits disagree due to difference in rounding off. 

The validity of integrals of Set IV relies on the fact that two versions of M§ agree com- 
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TABLE VI: Residual renormalization constants needed for the calculation of the mass-dependent 
contributions from Set IV diagrams. Notations are those of Eq. (|17p . 





0.000 721 65 (94) 




-0.000 018 910 (26) 


M (em) 


0.000 000 519 762 (21) 






aIjtS 6,P2 


0.000 705 (12) 


ALB^l 


-0.000 079 83 (10) 




0.000 255 64 (5) 




0.000 009 405 25 (83) 


A T ( em ) 


-0.000 000 779 612 (11) 




0.000 001 559 224 (19) 



AMgpj 5.374 (45) AM^ -0.628 832 • 



(me) 
2 

ALB^pl 1.476 3 (33) ALB { ™ p e) 2 -0.308 75 (32) 

ASm^pl 11.151 39 (32) ALB^l 1.885 733 (16) 

A4?P2 -1.641 436 (54) AS 2%P2 3 - 282 872 ( 107 ) 



M^p{ 1.094 258 •• 





0.038 01 (14) 




-0.001 641 9 (18) 


M (mt) 
1V1 2,P2 


0.000 078 067 4 (31) 








0.023 97 (29) 




-0.004 155 7 (45) 




0.015 483 (25) 


ALB^l 


0.000 831 107 (75) 


A T ( m *) 
^-^2*,P2 


-0.000 117 097 (15) 




0.000 234 (1) 



pletely with each other. As was noted in Sec. HI] we actually used only the second version of 
Mg to build integrals of tenth-order diagrams of Set IV, because of a technical problem in 
the first version. However, we are convinced that the integrals of Set IV are indeed bug-free. 

As is seen from (12T|) the contribution of Set IV to the muon g — 2 is sizable, which is not 
unexpected. This is because the order of magnitude of the contribution of the dominant 
(me) term can be readily estimated, noting that the leading ln(m M /m e ) term is determined 
by the charge renormalization procedure. This leads to 

4 10) [Set IV (me) ] ~ 4K 2 «( 8) [Group V] ~ -31.0 (24) 

where the factor 4 comes from the number of virtual photon lines of [Group V] into 
which a vacuum-polarization loop can be inserted, [Group V] ~ -2.179 (3) , and 
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TABLE VII: Contributions of diagrams Moi,p2 r • • , M24,p2 of Set IV to a M for (l\h) = (me). The 
multiplicity np is the number of vertex diagrams represented by the integral and is incorporated 
in the numerical value. The integrals M\2,P2, Mi6,P2; and Mxs,P2 are evaluated with quadruple 
precision. All other integrals are evaluated with double precision. 



Integral np Value (Error) Sampling per No. of 

including np iteration iterations 



01,P2 


28 


-0.369 


5 


(130) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


48 


M (me) 
u2,P2 


56 


-8.223 


2 


(308) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


140 


mM 

■* HQ DO 


28 


-3.986 


6 


(230) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


80 


mM 


56 


46.429 


2 


(587) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


156 


M (me) 

JVinc; DO 


56 


19.803 


9 


(105) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


48 


M (me) 

iflrif; DO 

M (me) 
07, P2 


56 


-11.614 


8 


(211) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


80 


56 


-0.558 


3 


(129) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


72 


08, P2 


56 


-48.877 


7 


(420) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


130 


M (me) 
JW 09,P2 


56 


4.817 


2 


(310) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


140 


M (me) 

ivilfl DO 


56 


18.291 


5 


(459) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


156 




28 


21.337 


7 


(299) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


100 




28 


-56.967 


6 


(174) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


20 


M (me) 

10, -PZ 


56 


-61.803 


8 


(142) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


48 


(me) 
JW 14,P2 


56 


21.147 


2 


(238) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


80 


M (me) 

JW 15,P2 


56 


7.639 


9 


(138) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


72 


M (me) 


56 


62.954 


8 


(414) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


35 


M (me) 


56 


62.823 


6 


(412) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


156 


M (me) 
JW 18,P2 


56 


-44.191 


1 


(207) 


1 


X 


10 8 , 


1 X 


10 9 


50, 


30 


(me) 
iw 19,P2 
M (me) 
iw 20,P2 


28 


-12.057 


1 


(14) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


40 


56 


9.281 


7 


(84) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


48 


M (me) 
JW 21,P2 


28 


4.359 





(12) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


40 


M (me) 
JW 22,P2 


56 


-2.934 


2 


(105) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


48 


M (me) 
iw 23,P2 


56 


-44.431 


4 


(185) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


72 


M (me) 
JW 24,P2 


56 


19.396 


5 


(175) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


72 



the enhancement factor 



ln(m M /m e 



3.6 . 



(25) 



The value ( 1241) may be regarded as a fair approximation to ( 1211 . 

By now we have evaluated the complete set of tenth-order diagrams containing vacuum- 
polarization subdiagrams {3, 2]]- 27]. (Note that the remaining Sets have no vacuum- 
polarization loop.) In particular its (me) contribution to the muon g — 2, namely all sets 



16 



TABLE VIII: Contributions of diagrams M 2 5 t p2,- ■ ■ , -/^47,P2 of Set IV to a M for (hh) 
multiplicity np is the number of vertex diagrams represented by the integral and is 
in the numerical value. All integrals are evaluated with double precision. 



Integral n f 



Value (Error) 
including uf 



Sampling per 
iteration 



= (me). The 
incorporated 



_ No. of 
iterations 





28 


-1.148 


1 


(73) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


40 


j/M 

iw 26,P2 




— 1 "3 ^79 


r. 
o 


(1 771 


i 


v 
a. 


10 8 


_L A 


10 9 


OU] 


80 

OU 


1V1 27,P2 


•56 


11.510 


8 


(246) 


1 


x 


10 8 , 


1 X 


10 9 


80, 


80 


M (me) 
IV1 2S,P2 


uU 


■^fi n 

OU.OIU 


n 

u 




1 

_L 


a 


10 8 


± A 


10 9 


OU, 


-LOU 


M (me) 


28 


36.421 


2 


(298) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


88 


30, P2 


28 


-43.675 


1 


(402) 


1 


x 


10 8 , 


1 X 


10 9 


80, 


156 


M (me) 
iw 31,P2 
M (me) 
1V1 32,P2 


28 


35.547 


o 


(22) 


1 


x 


10 8 , 


1 X 


10 9 


80 


40 


•56 


—30.768 


6 


(60) 


1 


x 


10 8 , 


1 X 


10 9 


80, 


40 


33, P2 

mH 

JW 34,P2 


28 


-14.328 


3 


fill 


1 


x 


10 8 , 


1 X 


10 9 


80, 


40 






Q 
o 


l UU 1 


1 

_L 


v 

A 


10 8 
±yj ? 


± A 


10 9 


so 

OU, 


40 

*±u 


M (me) 
JKi 35,P2 


•56 


—7.263 


8 


(65) 


1 


x 


10 8 , 


1 X 


10 9 


80, 


40 


M (me) 
JW 36,P2 


•56 


3.414 


7 


(871 


1 


x 


10 8 , 


1 X 


10 9 


80 


48 


M (me) 
iw 37,P2 


28 


7.820 


6 


(24) 


1 


x 


10 8 , 


1 X 


10 9 


80, 


40 


M (me) 
iw 38,P2 

iw 39,P2 


28 


-16.252 


5 


(81) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


48 


56 


-7.644 


5 


(78) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


40 


M (me) 
iw 40,P2 


56 


2.819 





(158) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


72 


j/W 

JW 41,P2 


28 


-25.777 


5 


(194) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


72 


M (me) 
iw 42,P2 


28 


26.504 





(291) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


88 


M (me) 
iw 43,P2 


28 


-29.852 





(100) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


48 


M (me) 
iw 44,P2 


56 


37.460 


8 


(173) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


72 


JW 45,P2 


28 


9.735 


1 


(173) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


72 


M (me) 
iw 46,P2 


28 


8.399 


9 


(213) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


80 


M (me) 

JW 47,P2 


28 


-22.410 


1 


(284) 


1 


X 


10 8 , 


1 X 


10 9 


80, 


104 



excluding light-by-light scattering loops, is given by 

4 10) [A11 sets excluding 1-1 loops] (me) ~ 48.88 (19). 



(26) 



This may be compared with the corresponding result A 



based on the renormalization group method 



35(]. 



'(10) 



32 obtained by an estimate 
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TABLE IX: Contributions of diagrams Moi ) p2,- • • , -^24,P2 of Set IV to a M for (hh) = (mi). The 
multiplicity uf is the number of vertex diagrams represented by the integral and is incorporated in 
the numerical value. The integral M\2 t p2 is evaluated with quadruple precision. All other integrals 
are evaluated with double precision. 



llltCgl Oil 


Tl J? 


Value (Error) 


Sampling per 


No. of 




including 


riF 


iteration 


iterations 


M (mt) 

lVlf\-\ DO 
Ul,r Z 


28 


-0.026 


89 


(15) 


1 


X 


10 7 


100 


M (mt) 

1VJ -C\0 P9 


56 


-0.002 


06 


(37) 


1 


X 


10 7 


100 


M (mt) 

lv± (X r K P9 


28 


0.051 


83 


(23) 


1 


X 


10 7 


100 


M (mt) 
04 P2 


56 


0.049 


99 


(55) 


1 


X 


10 7 


100 


M (mt) 

UO,r Z 


56 


-0.011 


20 


(29) 


1 


X 


10 7 


100 


M (mt) 

UD,r Z 

JW 07,P2 


56 


-0.060 


08 


(52) 


1 


X 


10 7 


100 


56 


-0.065 


38 


(24) 


1 


X 


10 7 


100 


M (mt) 

1V± (\Q P9 

M (mt) 

iUnn DO 


56 


-0.038 


51 


(55) 


1 


X 


10 7 


100 


56 


-0.026 


52 


(62) 


1 


X 


10 7 


100 


JW 10,P2 


56 


-0.050 


11 


(44) 


1 


X 


10 7 


100 


M (mt) 

±±,_/^Z 


28 


0.016 


61 


(16) 


1 


X 


10 7 


100 


JW 12,P2 


28 


-0.037 


85 


(2) 


1 


X 


10 7 


50 


M (mt) 


56 


-0.060 


34 


(30) 


1 


X 


10 7 


100 


{vat) 
1V1 IA,P2 


56 


0.000 


86 


(43) 


1 


X 


10 7 


100 


M (mt) 

JW 15,P2 

M (mt) 


56 


0.009 


38 


(23) 


1 


X 


10 7 


100 


56 


0.040 


26 


(49) 


1 


X 


10 7 


100 


M (mt) 


56 


0.101 


13 


(62) 


1 


X 


10 7 


100 


M (mt) 
IV1 1S,P2 

(mt) 
iw 19,P2 


56 


0.011 


93 


(28) 


1 


X 


10 7 


100 


28 


-0.063 


10 


(5) 


1 


X 


10 7 


100 


M (mt) 

JW 20,P2 

M (mt) 

m 2\,P2 

M (mt) 

IVI 22,P2 

M (mt) 

IV1 23,P2 

M (mt) 

M 24,P2 


56 


0.029 


21 


(37) 


1 


X 


10 7 


100 


28 


0.007 


11 


(4) 


1 


X 


10 7 


100 


56 


-0.034 


71 


(33) 


1 


X 


10 7 


100 


56 


-0.008 


27 


(68) 


1 


X 


10 7 


100 


56 


0.022 


23 


(39) 


1 


X 


10 7 


100 
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TABLE X: Contributions of diagrams M25 ) p2,. • • , -^47,P2 of Set IV to for (hfo) = (mt). The 
multiplicity np is the number of vertex diagrams represented by the integral and is incorporated 
in the numerical value. All integrals are evaluated with double precision. 



Integral 


np 


Value (Error) 


Sampling per 


INo. OI 




including 




iteration 


iterations 


lyl 2h,P2 


28 


-0 039 19 


(15) 


1 X 


10 7 


100 


M (mt) 
1V1 26,P2 


28 


— D DID 15 


(17) 


1 X 


10 7 


100 


1V1 27,P2 


56 


060 33 


(53) 


1 X 


10 7 


100 


Jla 28,P2 
1V1 29,P2 


56 


006 03 


(36) 


1 X 


10 7 


100 


28 


037 86 


(34) 


1 X 


10 7 


100 


30, P2 


28 


-0.005 19 


(24) 


1 X 


10 7 


100 


31, P2 


28 


0.056 43 


(12) 


1 X 


10 7 


100 


32, P2 


56 


-0.096 26 


(30) 


1 X 


10 7 


100 


33, P2 
M (mt) 
JW 34,P2 


28 


-0 030 72 




1 X 


10 7 


100 


56 


-0 060 57 


(25) 


1 X 


10 7 


100 


M (mt) 


56 


-0 025 98 


(31) 


1 X 


10 7 


100 


Mir! 

JW 36,P2 


56 


050 76 


(32) 


1 X 


10 7 


100 


M (mt) 


28 


019 70 


(10) 


1 X 


10 7 


100 


M^l 
38, P2 


28 


021 52 


(14) 


1 X 


10 7 


100 


M (mt) 
JW 39,P2 


56 


-0.034 44 


(21) 


1 X 


10 7 


100 


M (mt) 
JW 40,P2 


56 


-0.018 44 


(21) 


1 X 


10 7 


100 


M (mt) 

M (mt) 

1V1 A2,P2 

M (mt) 


28 


-0.033 08 


(28) 


1 X 


10 7 


100 


28 


0.006 74 


(25) 


1 X 


10 7 


100 


28 


-0.016 64 


(36) 


1 X 


10 7 


100 


M (mt) 

M (mt) 
1V1 45,P2 


56 


0.039 39 


(35) 


1 X 


10 7 


100 


28 


0.002 23 


(14) 


1 X 


10 7 


100 


M (mt) 
1VI 46,P2 


28 


-0.028 76 


(34) 


1 X 


10 7 


100 


M (mt) 
IVI 47,P2 


28 


-0.004 02 


(21) 


1 X 


10 7 


100 
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Appendix A: Summing up Residual Renormalization Terms 

The purpose of this Appendix is to obtain the sum of residual renormalization terms of 
the Set IV. Since diagrams of Set IV have exact correspondence with the diagrams of Group 
V of the eighth- order g — 2, however, it is simpler to consider the residual renormalization 
of the diagrams of Group V, from which the residual renormalization of the Set IV can be 
readily derived. 

In our approach integrals of individual diagrams must be made convergent before they are 
integrated numerically. This is achieved by constructing terms which subtract UV-divergent 
parts by A-operation and IR-divergent parts by //^-subtractions. Since this scheme is 
different from the standard on-shell renormalization, it is necessary to make an adjustment, 
called residual renormalization. Residual renormalization terms of individual diagrams must 
then be summed up over all diagrams involved. 

As the order of perturbation increases the total number of terms contributing to the resid- 
ual renormalization increases rapidly so that it will become harder and harder to manage. 
Fortunately the sum of all residual terms can be expressed concisely in terms of magnetic 
moments and finite parts of renormalization constants of lower orders j^J , and the sum has 
a structure closely related to that of the standard on-shell renormalization. This enables 
us to confirm the validity of the sum of residual renormalization terms starting from the 
expression of the standard renormalization. 

To see this relation clearly it is useful to treat UV- divergence and IR-divergence sepa- 
rately. We present the logic of our approach for the fourth-, sixth-, and eighth-order cases, 
in that order. We deal here only with Ward-Takahashi(WT)-summed diagrams of g-type, 
namely diagrams without closed lepton loops. Thus M 2n and ami n — 1> 2, • • • , refer to 
unrenormalized and renormalized amplitudes of such diagrams, respectively. 

Our discussion here follows the scheme incorporated in the automatic code generator 
GENCO-deA, which is applicable to any value of the order A. 
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1. fourth-order case 



The standard renormalization of the fourth-order magnetic moment a 4 can be expressed 
in the form 

a 4 = M 4 - 2L 2 M 2 - B 2 M 2 - Sm 2 M 2 * , (Al) 

where M 2 is the second-order magnetic moment, M 2 * is obtained from M 2 by inserting a 
two-point vertex in the lepton line of M 2 , and M 4 is the sum of unrenormalized WT-summed 
amplitudes M 4ct and M^: 

M 4 = M 4a + M 4b , (A2) 

where 4a and 46 refer to the fourth-order diagrams in which two virtual photons are crossed 
and uncrossed, respectively. The coefficients of renormalization constants L 2 and B 2 in 
Eq. (lAip reflect the fact that M 4a is obtained by inserting a second-order vertex diagram in 
two vertices of M 2 and is obtained by inserting a second-order self-energy diagram in 
the electron line of M 2 . 

a. Separation of UV divergences by the K-operation 

M 4 has no overall UV divergence. However, it has UV divergences coming from subdia- 
grams. Applying i^-operation on these divergences we obtain 

M 4 = B^ V M 2 + 5m^ y M 2 * + 2L% V M 2 + M 4 R , (A3) 

where the superscript R in Mf means that all subdiagram UV divergences are removed 
from M 4 . and B^ Y are the UV-divergent parts separated out from L 2 and B 2 by the 
.ff-operation and and Bf^ are UV-finite (but IR-divergent) remainders: 



L 2 


— + 


B 2 


= B^ + B* 


5m 2 


= 5mf. 



(A4) 



5m\^ = is the specific feature of the i^-operation for the second-order self-energy diagram. 
Substituting Eqs. (|A3]) and (TX4]) in Eq. flAT]) we obtain 

a 4 = Mf - M 2 (2Lf + Bf) , ( A5) 
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Note that the coefficients of L R and _B R in Eq. (1A5j) inherit the coefficients of L2 and B2 in 
Eq. fl5H). 



6. Separation of IR divergences by the I /R- subtraction 

The second-order mass renormalization is completely carried out and no remainder is 
left in the i^-operation. The .R-subtraction, then, is not applied by GENCODE./V in the 
case of the fourth order. IR divergence is caused by a photon spanning over a self-energy- 
like subdiagram which actually represents a lower-order magnetic moment. This magnetic 
moment can be effectively represented by a three-point vertex between one photon and two 
electrons. Thus, the UV-fmite term M R must have an IR singular structure which is similar 
to that of the vertex renormalization constant L R : 

Mf = M 2 Lf + AM 4 , (A6) 

where M2 comes from the second-order self-energy subdiagram of and L R appears by 
replacing the M2 self-energy subdiagram by a point vertex. 

The IR-divergence is also found in the vertex and wave-function renormalization con- 
stants. The WT-identity 

L 2 + B 2 = (A7) 

guarantees that L 2 and B 2 have the same, but opposite in sign, IR singularity. This enables 
us to separate the IR-singular and finite terms of L R and _B R as follows: 

L R = I 2 + AL 2 , 

Bf = -I 2 + AB 2 . (A8) 
where I2 is IR-singular but its finite part is undetermined. The finite terms AL2 and AB2 



depend on how we define J2. For instance, in Ref. j32j, the /-operation was defined so that 
I2 = £ R = ln(A/m) + 5/4, where A is the photon mass. The sum L R + Bf, however, does 
not depend on the definition of 12- We find that 

ALB 2 = L^ + Bf = AL 2 + AB 2 = ^. (A9) 

In other words, the finite quantity ALB2 is determined by how we extract UV divergence 
by the i^-operation from each of L2 and B2. 

L^ + B^ = - 3 -. (A10) 
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Substituting Eqs. (1A6j) and ( 1A8j) in Eq. ( 1A5j) . one can express 04 denned by the standard 
renormalization as a sum of finite terms only: 

a 4 = AM 4 — M 2 ALB 2 . (All) 

2. sixth-order case 

The sixth-order magnetic moment has contributions from ten diagrams, each of which 
represents the sum of five vertex diagrams transformed with the help of the WT-identity. In 
the standard renormalization it can be written in terms of unrenormalized amplitudes Mq, 
M4, etc., and various renormalization constants as 

a 6 = M 6 

- M 4 (3B 2 + 4L 2 ) - M 4 *5m 2 

- M 2 (B 4 + 2L 4 ) - M 2 *Sm 4 

+ M 2 {2(B 2 ) 2 + 8B 2 L 2 + 7(L 2 ) 2 } + M 2 *(3B 2 + AL 2 )5m 2 
+ M 2 {B 2 * +AL 2 *)5m 2 

+ M 2 *Sm 2 5m 2 * + M 2 »»(5m 2 ) 2 , (A12) 

where M 4 is defined by Eq. (IA2[) . M 2 »» is obtained from M 2 by inserting two two-point 
vertices in the lepton line of M 2 , and 

H 

M 6 = y^jMei, rji = l except that r] D = r] G = 2, 

i=A 

3 

M 4 * = ^(M 4a(i * } +M 46( ^)), 

i=l 

B 4 = B 4a + B^, 

3 

L4 = ^J(jy4aQ) + L/tb(i)), 
i=l 

5m 4 = <5m 4a + <5m 4b . (A13) 

Mj a (i*) is obtained from M 4a (which contains three lepton lines 1, 2, 3) by inserting a two- 
point vertex in the lepton line i of M 4a , and L 4a (j) is the vertex renormalization constant of 
the diagram in which an external vertex is inserted in the lepton line i of the diagram 4a. 
Similar notation is applied for the diagrams built from 4b. 
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The coefficient of M 4 in Eq. (1A12j) can be readily understood noting that the fourth- 
order self-energy diagrams M 4a and M 4 & have three fermion lines into which second-order 
self-energy can be inserted and four vertices into which second-order vertex can be inserted. 
Similarly, there are one fermion line and two vertices in the second-order self-energy diagram 
M 2 into which we can insert a B4 or a L 4 , which leads to — M 2 (5 4 + 2L 4 ). The term 
M 2 {2(5 2 ) 2 + 8B2L2 + 7(L2) 2 } comes from two ways of inserting B2 in M2 (one disjoint and 



one nested relations of two _B 2 's 



19|). eight ways of inserting one L2 and one B2 in M2 



(two disjoint, two overlapping, and four nested relations of L 2 and B 2 ), and seven ways of 
inserting two L 2 in M 2 (one disjoint, four overlapping, and two nested relations of two L 2 's). 
There is only one way to insert <5m 4 in M 2 and 5m 2 in B 2 of M 2 B 2 . There are three ways 
to insert 5m 2 in M 4 , but the coefficient three is included in the definition of M 4 *. There are 
four ways to insert 5m 2 in L 2 of M 2 L 2 . The coefficients of other terms can be understood 
in a similar fashion. 



to 



a. Separation of UV divergences by the K-operation 

Analysis of the UV divergence structure of M 6 , L 4 , £? 4 , and 5m 4 by the if-operation leads 

M 6 = M 6 R 

+ M 4 (35 2 uv + 4L 2 JV ) + M 4 *5m 2 

+ M 2 (B™ + 2L 4 JV ) + M 2 *5m 4 JV 

- M 2 {(S 2 UV ) 2 + BffB™ + AB^Lf + AB^L^ + 7(L^ y ) 2 } 

- M 2 *(2£ 2 JV + 4L 2 JV )5m 2 

- M 2 *B™ 6m ™ 

- M 2 *5m 2 7Sm 2 

- M 2 **(6m 2 ) 2 , (A14) 



where 



M 6 R = J2^ iM ^ = 1 except that rj D = Vg = 2, 



i=A 



(A15) 
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is the UV-finite part of Mq. UV-divergent parts of L 4 , B4, and £m 4 are separated as follows: 
L 4 = Z^v + 3L^ Y Lf + 25 2 uv Lj + 25m 2 L 2 * + 

<5m 4 = SmY* + 5m 2 5m% + B^ y 5m5 + 5m*. (A16) 



B 4 = BY V + 2L^B R + B™B* + 5m 2 B 2 « + B? 



Mf is denned in Eq. (jA"3]) . and L^ v and are denned in Eq. flAij) . 

Substituting Eqs. flATil) . flA3l) . flAT6l) . and (jA4|) in Eq. flAT2l in this order, we obtain a 6 
expressed by UV-finite quantities only: 

a 6 = M 6 R 

- M 4 R (35* + 

- M 2 {Bf + 2L*)-M 2 *5m* 

+ M 2 {2(B*) 2 + SBfLf + 7(L*) 2 } . (A17) 

Note that this equation has exactly the same structure as Eq. (1A12I) . although it looks 
simpler because 5m* = in the i^-operation. This is what one would expect since, in 
Eq. (1A12[) . all UV-divergent quantities must cancel out, leaving only UV-finite pieces with 



their original numerical coefficients. 



b. Separation of IR divergences by the I /R- subtraction 



Since Eq. (1A17|) has no linearly IR divergent term caused by the self-mass term, there is no 



need to invoke the i?-subtraction. We, however, retain the i?-subtraction that is incorporated 
in GENCODEiV. Quantities obtained above can be expressed as the sum of logarithmically 
IR-divergent pieces defined by the /-subtraction and finite remainders together with the 
residual mass-renormalization term defined by the i?-subtraction: 

Mf = L R M 2 - (L K ) 2 M 2 + LfMf + 5mfM 2 , + AM 6 , 
L* = / 4 + (L*) 2 + AL 4 

Bf = -I 4 + LfB* + AB A , (A18) 

where IR-divergent terms are contained in L*, B* and J 4 term. The WT-identity guarantees 
that L 4 and B4 have the same overall IR-divergence which we call 7 4 . In the previous work 
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32J the 7 4 is chosen as the sum of non-contraction terms I 4a (i) of the vertex renormalization 
constants L ia ^y. 

h = ha(l) + ha(2) + ha(3) + hb(l) + ^46(2) + ^46(3) ■ (A19) 

The finite quantities AL 4 and A5 4 depend on how J 4 is defined. But the sum of + Bf is 
independent from the definition of J 4 . Therefore, we introduce the finite quantity ALB 4 by 

ALB4 = L 4 + £ 4 - LfAW 2 = AL 4 + A5 4 • (A20) 

Note that the value of ALB 4 is unambiguously determined by our choice of L^ v and 
in the i^-operation and by the WT-identity L 4 + B<± = 0. 

Substituting Eqs. (IA18p . (IA20p . dA6|) . and (JA9]) in Eq. flA17j> in this order, we obtain a 6 



of standard renormalization as the sum of finite terms only 

a 6 = AM 6 - 3AM 4 ALB 2 

+ M 2 {-ALB 4 + 2(ALS 2 ) 2 }. (A21) 

3. eighth-order case 

The eighth-order magnetic moment a 8 has contributions from 74 WT-summed diagrams. 
In the standard renormalization the renormalized moment a 8 can be written in terms of 



26 



unrenormalized amplitudes M 8 , M 6 , M 4 , etc., and various renormalization constants as 

a 8 = M 8 

— M 6 (55 2 + 6L 2 ) 

— M 6 *Sm 2 

+ M 4 {-3B 4 - 4L 4 + 9(B 2 ) 2 + 2QB 2 L 2 + 18(L 2 ) 2 + 5m 2 (3B 2 * + 8L 2 *)} 
+ M 4 *{5m 2 {hB 2 + 6L 2 ) + Sm 2 Sm 2 * — 5m 4 } 
+ M 4 **(5m 2 ) 2 

+ M 2 {-B 6 - 2L 6 + \2L 4 B 2 + 18L 4 L 2 + 654^2 + WB 4 L 2 

-545 2 (L 2 ) 2 - 30(5 2 ) 2 L 2 - 5(5 2 ) 3 - 30(L 2 ) 3 } 
+ M 2 5m 4 (B 2 * + 4L 2 .) 

+ M 2 6m 2 (B 4 , + 2L 4 * - 6B 2 B 2 * - 24B 2 L 2 * - 10B 2 *L 2 - 36L 2 L 2 ») 

— M 2 5m 2 5m 2 *(B 2 * + 4L 2 * ) 

— M 2 (5m 2 f(B 2 „ +4L 2 » t + 2L 2 , t *) 

+ M 2 «5m 2 {3B 4 + 4L 4 + 5m 4 * - 26B 2 L 2 - 9(5 2 ) 2 - 18(L 2 ) 2 } 

— M 2 *5rriQ 

+ M 2 *5m 4 (5B 2 + 6L 2 + 5m 2 *) 

— M 2 *(5m 2 ) 2 (3B 2 * + 8L 2 * + 5m 2 **) 

— M 2 *8m 2 5m 2 *{hB 2 + 6L 2 ) 

— M 2 *5m 2 (5m 2 *) 2 

+ M 2 **5m 2 {25m 4 — 5m 2 (5B 2 + 6L 2 + 2<5m 2 .)} 

— M 2 ,,,{5m 2 f. (A22) 
M 8 is defined by 

47 

M 8 = VaM a , (A23) 

Q = 01 

where r\ a = 1 for time-reversal-symmetric diagrams and r] a = 2 for others. 
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a. Separation of UV divergences by the K-operation 



The UV divergence structure of M 8 is given by 

M 8 = M 8 R 

+ M 6 (5B™ + 6L™) 
+ M 6 .<5m 2 

+ M 4 {35 4 uv + 4L™ - 3B™B™ - 6(5 2 uv ) 2 - 18B™L™ - 8B™L™ _ 18(L 2 JV ) 2 } 
+ M 4 .(<5mY V - B^Sm^ - 45m 2 5 2 uv - 6Sm 2 L™ - Sm 2 Sm^7) 
— M 4 ..(5m 2 ) 2 

+ M 2 {B™ + 2L^ V - 2B™B™ - 6S 4 UV LU V - B^B™ - 4S 4 UV L 2 J , V 

-4L? V B 2 UV - lSL^™ + GB^L^B^ + 36B™L™L™ + 18B™(L™) 2 
-B™B™ - 2B™L™ + 45 2 UV J B 2 U V L 2 J , V + £ 2 uv (i^ v ) 2 + 6( J B 2 UV ) 2 L 2 JV 
+2{B™) 2 L™ + (B™) 2 B™ + 2(B™) 2 B™ 
+8(B™) 2 L™ + (£ 2 UV ) 3 + 30(L 2 JV ) 3 } 

+ M 2 *{6m™ -B™5m™ + (B™) 2 5m™} 

+ M 2 »5m™{-2B™ -6L™ -6m™} 

+ M 2 *5m™{-B™ + 6B™L™ + 2(B uv )2 + ^uv^uv + B ^ S m™} 
+ M 2 *5m 2 {-2B™ - AL™ - tfm™ + 18(L 2 V ) 2 

+B™(12L V 2 Y + 8L™ + 2B™ + 35 2 uv + 25m™) 

+5m™(2B™ + 6L™ + 5m™)} 
+ M 2 **5m 2 {-25m™ + (35m 2 + 25m™)B™ + Q5m 2 L 2 ]V + 25m 2 5m™} 
+ M 2 ***(8m 2 f . (A24) 

The quantities with a prime, L 2 >, called derivative amplitudes, arises from a fourth-order 
diagram that contains a self-energy subdiagram. This self-energy subdiagram supplies the 
inverse of the fermion propagator times the wave function renormalization constant and can- 
cels out one of the adjacent fermion propagators and yields another renormalization constant 
L 2 . In the expression L 2 >, the inverse fermion propagator and the adjacent fermion propa- 
gator are left in the numerator and denominator, respectively, of the Feynman parametric 
expression of the amplitude. Thus the whole renormalization constant L 2 > is analytically 
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identical with L 2 . But, the separation of UV divergence by the if-operation works differently 
in two cases, so that L^, v is different from by a finite amount. 

A similar consideration applies to higher order quantities. Take, for instance, L 4 which 
consists of four fermion lines. There are four ways to insert a self-energy subdiagram to 
L 4 . Since we define L 4 > as the sum of all derivative amplitudes of L 4 , we have L 4 > = 4L 4 . 
Similarly, B 4 > = 3B 4 . 

The second-order derivative amplitude, such as L 2 /, however, does not include its sym- 
metric factor in our definition. Thus, L 2 > = L 2 . 

We also need the UV divergence structures of the renormalization terms B 6 , L 6 , and 5m 6 : 




+ B 4 *Sm 2 



+ Bf{2L v ^ - 45 2 UV LU V - 7(LU V ) 2 } 




- bUbT? 

+ B^,(-25m 2 B^ v ) 



+ B A ,{B^) 



+ B 2 *(5m){ v - B2 5m 2 , - 45m 2 L^ 
- B 2 **(5m 2 ) 2 , 



Sm 2 Sm 2 , ) 



(A25) 
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and 



t _ rUV , rR 

-^6 — -^6 -^6 

+ L A *(8m 2 ) 

+ Lf{3L^ Y - QB^iy/ - 12(LU V ) 2 } 

+ L£(2£ 4 uv - lOS^L™ - 2£ 2 UV £# V ) 

rR / p>UV\2 

+ Z^L^) 

+ M£ 2 UV ) 

+ L R ,(-25m 2 5 2 uv ) 

+ L 2 *{28m^ Y - 2B% v 6rr$ / - 108m 2 L 2 JY - 28m 2 8m^) 

- L 2 **(8m 2 ) 2 , (A26) 



- 5m R (B 2 UV ) 2 

+ 6rr$(BY v -B™B$ v ) 

+ 5mf,{5m^ y - B^ Y 5m^ v - 8m 2 8m^ 

+ 8m%,,(-28m 2 B$ v ) 

+ 8rn^{AL\ JV ) 

+ 5m A *5m 2 



+ 5m A ,B VY 



2 

- 5m 2 „(5m 2 f, (A27) 

where the quantity A is defined hy A = A — A VY . The difference between A and A R is 
that the former contains UV divergent terms arising from subdiagrams, while the latter is 
completely free from these sub-UV divergences. For instance, we have 

pT r r uv 

= 5 4 R + <5m 2 5 2 * + B^B^, + 2L 2 f v 5 R (A28) 

= L R + 25m 2 L 2 * + 2B^ v Lf, + 2L 2 JV L R (A29) 
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and so on. 

We also need the UV divergence structure of M 4 «, which is the amplitude of the fourth- 
order magnetic moment with a two-point vertex insertion: 

M 4 * = M R + 2Lf / M 2 * + 2(5m 2 M 2 ,, + B% v M 2 .) + 5m%? ' M v . (A30) 

Substituting the UV structures of the eighth order Eq. fl A24j) . the sixth-order quantities 
Eqs. (TATijl . (TA25]) . (TA26]) and flA"27]) . those of the fourth order Eqs. (TM]) . (lATBj) . and (TMD]), 
and those of the second order (IA4[) in this sequence in Eq. flA22p . we obtain the UV- finite 
expression of the magnetic moment a§: 

a 8 = M R 

+ M*(-5Bf - 6L R ) 

+ M 4 R {-4L R - 3Bf + 26L R 5 2 R + 18(L R ) 2 + 9(£ 2 R ) 2 } 
- M R 5m R 



+ M 2 {-B 6 R - 2L R + (fl 2 . + 4L 2 ,)Sm 



+6BfBf + lOBfLf + 125 R L R + 18L R L R 
-30L R (5 R ) 2 - 54(L R ) 2 J B R - 30(L R ) 3 - 5(5 R ) 3 } 
+ M 2 »{-<5m R + <5m R (5m R + QLf + 5Bf)} . (A31) 

Again Eq. (IA31j) has exactly the same structure as Eq. (IA22j) except that #m R = 0. 



6. I /R- subtraction 

In order to handle the numerical calculation on a computer, we need to separate the 
IR divergent terms from M R . Paying attention to the outermost photon spanning over a 
self-energy subdiagram, we obtain the IR structure of M% as follows: 

Mf = AM 8 
+ M 6 R L R 

+ Mf{Lf - (L R ) 2 } 

+ M 2 {L R - 2L R L R + (L R ) 3 - 25m R L R } 
+ M R 5m R 

+ M 2 ,(5m R - 5m R 5m R - 5m R L R ) . (A32) 
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Eq. ( 1A32I) has a term Mf,5m^, where M$ is linearly IR-divergent, which arises from the 
diagrams M±e and Mig. This term compensates the same IR-divergence found in a 8 of 
Eq. (IA31j) whose origin is the mass-renormalization term M^Sm^ associated with the dia- 
grams Mie and Mi 8 . This IR divergence in Mf- can thus be removed from M\& and Mig 
by the i?-subtraction which acts on a fourth-order self-energy subdiagram of M 16 (M 18 ) and 
complements the mass-renormalization constant Sm^Sm^) . 

Another linear IR divergent term in Eq. f)A32|) is 2M 2 L^,5m^ , where L 2 . is linearly 
divergent. This IR divergence is again found in the diagrams of Mi$ and Mi 8 . In the IR-limit 
of the outermost photon loop, a possible configuration of Mi 6 (Mi 8 ) is that the second-order 
self-energy subdiagram of M\q(Mis) supplies the second-order anomalous magnetic moment 
M 2 and the fourth-order self-energy subdiagram behaves as Srn^uy The IR behavior of 
the residual diagram including the outermost photon line resembles the second-order vertex 
diagram with a two-point vertex insertion L\\ : 

L 2 * = AL 2 * + L\\ , (A33) 

where AL 2 . = —3/4 is the one contraction term of L 2 » and the IR divergent L\\ is the 
non-contraction term of L 2 *. 

This IR divergence in Mf- of Eq. flA32j) compensates the IR divergence in 2M 2 L 2 *5m4 
of the renormalized magnetic moment a 8 of Eq. f ]A31j) . The origin of the +4L 2 *#m4M 2 in 
Eq. flA3ip is the renormalization terms associated with the diagrams M 08 , M w , M 41 , and 
M46. Two of four L 2 * terms are exactly canceled by _B 2 * terms because of the WT-identity 
2L 2 . + B 2 * = 0. The remaining two L 2 . terms will cancel the IR-divergence arising from the 
diagrams M 16 and M 18 in Mf". 

The last of the linearly IR divergent terms of of Eq. (1A32j) is M 2 L^, which also 



comes from M 16 and M 18 . In this case, the second-order self-energy subdiagram supplies 
a second-order anomalous magnetic moment M 2 and the rest of the residual diagrams are 
pushed in the IR limit. From Mis(Mig), it gives rise to the similar IR behavior of the six- 
order vertex renormalization constant L 6 6( 1 )(L 6c ( 1 )). This IR divergence will be canceled in 
a 8 of Eq. (1A31I) by the M 2 L^ term which comes from the renormalization constants for the 
diagram M 08 (M 10 ). 

To see the cancellation of remaining logarithmic IR divergence in a 8 , we need the IR 
structures of the renormalization constants Lf and Bf. Resorting to the WT-identity again, 
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we can define the finite quantity ALBq as follows: 



ALB, = L% + Bf 

- {+I 6 + 2LfLf-(Lf) 3 + 25m*L 2 «} 

- {-h + L^Bf + L*Bf - (Lf) 2 B* + 6mfB 2 * } 



(A34) 



where Iq is the overall IR divergent term of L 6 and Bq. The WT-identity guarantees that 
ALBq is independent of the choice of Iq. Note that ALBq ee AL 6 + ABq + AL^AB 2 + 
A<5m 4 £> 2 » [7] , where the quantities in the right-hand side are defined in Ref. [321 ] . 



c. Finite expression 

Separating the UV-finite quantities in a 8 of Eq. (1A31j) into the IR-singular parts and the 
finite parts as given in Eqs. (jA"32il . flAT8|) . flA6|) . flA34l) . flA20|) . and (jA"9l) . we obtain the 
expression a§ in terms of the finite quantities only: 

a 8 = AM 8 

+ AM 6 (-5ALB 2 ) 

+ AM 4 {-3ALB 4 + 9(ALB 2 ) 2 } 

+ M 2 {-ALB 6 + 6ALB4ALB2 - 5(ALS 2 ) 3 } 

+ 2M 2 AL 2 .A5m 4 . (A35) 
Since ALB 4 = AL 4 + AS 4 , 2AL 2 . = -AB 2 *, and ALB 2 = A52, this is equivalent to 



Eq. (76) of Ref. [5|, which was obtained from the direct sum of all subtraction terms. Note 
that the last term of Eq. (IA35I) remains unsubtracted regardless of the i?-subtraction, which 
is the residual mass-renormalization. This is because we use only the non-contraction term 

as the IR-subtraction term, leaving the finite part of AL 2 . untouched. 

The definition of the finite term AL 2 . does depend on how to separate IR part from L 2 * . 
To avoid such arbitrariness, we stick to the same /-subtraction rule of IR separation which 
is used for vertex renormalization constants. Namely, the IR-singularity is confined in L^, 
which is defined by the rule 

L n = L n — highest contraction term of L n 

ee L n — UV sub-divergence term determined by if-operation on L n (A36) 
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and this is used as an IR-subtraction term. This determines AL2* = —3/4 unambigu- 
ously. The iT-operation does not pick up this AL 2 * term from a corresponding subdiagram, 
since L 2 * is UV-finite. So, no rule exists in the automation code GENCODEiV that allows us 
to subtract the finite term AL 2 * of a renormalization constant. 

The residual renormalization scheme for the Set IV contribution [Set IV^ 1 ' 2 - 1 ] can be 
readily obtained from Eq. (1A35I) by insertion of a closed loop of the lepton / 2 in the internal 
photon lines of a 8 . This leads to Eq. (JT7J) given in Sec. IIHI 
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